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ABSTRACT 


Presents  surface  area  equations  for  ponderosa  pine  ("yellow" 
and  "blackjack"  pine  both  separately  and  in  combination) , 
Douglas-fir ,  and  a  combination  of  Mexican  white  pine  and  whiter- 
fir.    The  basic  data  came  from  three  National  Forests  in  Arizona 
and  New  Mexico,  and  were  analyzed  using  weighted  least  squares 
regression  techniques.    Also  presented  are  equations  for  ^ 
predicting  surface  area  of  forked  trees.  „ 
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The  total  stem  surface  area  of  a  tree  is  the  area  of  the  bole  surface  inside  bark, 
excluding  the  end  areas,  from  the  ground  to  the  tip  of  the  tree.     It  is  equivalent  to 
the  area  of  a  sheet  of  paper  wrapped  around  a  form  which  approximates  a  paraboloid  or 
cone.    The  use  of  surface  area,  when  combined  with  cubic-foot  volume  and  length  (tree 
height) ,  to  predict  product  recovery  from  trees  and  logs  has  been  described  by  Grosen- 
baugh  (1954),  Davis  and  others  (1962),  and  Bruce  (1970).     Lexen  (1943)  and  Husch  (1963) 
suggested  total  stem  surface  area  per  acre  as  a  measure  of  growing  stock  for  estimating 
stand  growth.     Stand  surface  area  inside  bark  is  equivalent  to  area  of  the  lateral 
meristem,  where  volume  growth  occurs. 

The  lack  of  surface  area  equations  which  could  be  used  in  product  recovery  studies 
in  the  Southwest  stimulated  the  research  reported  here.    This  paper  describes  methods 
used  to  develop  surface  area  equations  as  a  function  of  total  tree  height  and  d.b.h. 
for  four  species  (one  with  two  vigor  types)  or  species  groups. 
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SOURCE  AND  NATURE  OF  DATA 


The  data  were  collected  by  field  crews  from  the  Division  of  Timber  Management, 
Southwestern  Region,  USDA  Forest  Service,  located  in  Albuquerque,  New  Mexico. 

Species  and  Forests  were: 


Douglas-fir  (Pseudotsuga  menziesii  var.  glauaa  iMirb.]  Franco)  from  the 
Lincoln  and  Tonto  National  Forests. 

Mexican  white  pine  (Firms  flexilis  var.  veflexa  Engelm.)  from  the  Lincoln 
National  Forest. 

White  fir  (Abies  concolor  [Gord.  and  Glend.])  from  the  Lincoln  National 
Forest . 

"Blackjack  pine"  (Pinus  ponderosa  Laws.)  from  the  Lincoln,  Coconino,  and 
Tonto  National  Forests. 

"Yellow  pine"  (Pinus  ponderosa  Laws.)  from  the  Lincoln,  Coconino,  and 
Tonto  National  Forests. 


Foresters  in  the  Southwest  separate  "blackjack"  and  "yellow"  pines  on  the  basis 
of  growth  vigor.    The  distinction  is  meaningless  botanical ly  but,  because  of  the 
differences  among  trees  and  stands  of  low  and  high  growth  vigor,  it  is  useful  in  forest 
management . 

This  study  was  coordinated  with  the  continuous  forest  inventory  (CFI)  program  in 
the  Southwestern  Region.    A  subsample  of  the  10-point  clusters  measured  in  the  CFI 
program  was  randomly  drawn  for  each  Forest.     On  the  Coconino  and  Tonto  National  Forests, 
all  of  the  trees  located  on  the  10-point  cluster  were  felled  and  measured  using  pro- 
cedures similar  to  those  described  by  Stage  and  others  (1968) .     Because  of  a  lack  of 
funds,  and  for  safety  reasons,  not  all  of  the  trees  were  felled  on  the  Lincoln  National 
Forest.     Instead,  trees  were  selected  to  insure  representation  from  all  of  the  diameter 
classes  in  the  cluster. 


Sample  tree  surface  areas  were  calculated  with  a  modified  version  of  the  NETVSL 
computer  program  (Stage  and  others  1968).     Before  the  output  from  NETVSL  was  examined, 
expected  model  forms  were  developed.    The  expected  model  forms  were  used  to  recognize 
and  predict  the  relationship  between  surface  area  and  the  independent  variables  of 
total  tree  height  and  d.b.h. 

Three  simple  geometric  forms  were  considered.    A  tree  stem,  or  portions  thereof, 
can  be  approximated  by  a  cylinder,  cone,  or  paraboloid.    The  surface  area  for  each 
geometric  solid  can  be  expressed  by  the  following  formulas  (Husch  1963) : 


EXPECTED  MODEL  FORMS 


FORM 


FORMULA 


Cylinder 
Cone 

Paraboloid 


SA  =  27Trh 

SA  =  Tir(r2  +  h2)1/2 

SA  =  [(irr)/(6h2)]  [  (4h2  +  r2)3/2  -  r3] 


where  SA  =  surface  area 


r  =  radius  of  the  base,  in  feet 
h  =  height  of  the  form,  in  feet 
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Although  the  surface  area  formulas  for  the  cone  and  the  paraboloid  appear  complex, 
simplified  formulas  can  be  derived  without  great  loss  of  accuracy.     For  tree  stems,  the 
radius  (when  expressed  in  feet)  is  small  compared  to  height.    Therefore,  when  the  height 
and  radius  are  squared  (or  cubed)  and  added  (or  subtracted)  together,  the  radius  con- 
tributes little  to  the  result.    Hence,  the  two  formulas  were  simplified: 

FORMULA  FOR  CONE  FORMULA  FOR  PARABOLOID 

SA  =  Trr(r2  +  h2)1/2  SA  =  [  (irr) /6h2) ]  [  (4h2  +  r2)3/2  -  r3] 

SA  *  TrrQi2)1/2  SA  »  [Ur)/6h2)][(4h2)3/2] 

SA  =  ?rrh  SA  =  4/3  m:h 

From  this  simplification,  it  was  apparent  that  the  major  difference  among  the  surface 
area  formulas  was  the  constant  of  2tt  for  cylinders,  4/3tt  for  paraboloids,  and  it  for 
cones.    The  expected  model,  therefore,  is  the  linear  relationship: 


Surface  Area  =  bQ  +  bj (DH) 

where  D  =  diameter  at  breast  height,  in  inches 

H  =  total  tree  height,  in  feet 
b    =  an  intercept  which  should  range  from  0  to  4.0  square  feet.  An 
intercept  can  be  expected  because  trees  4.5  feet  tall  or  less 
have  DH  =  0. 

bj  =  a  constant,  which  should  range  from  — (24)    to  7A  ^or  ^rom 

0.1745  to  0.1309)  because  trees  are  between  a  cone  and  a 
paraboloid  in  shape. 

ANALYSIS  OF  UNFORKED  TREE  DATA 


Surface  area  of  unforked  trees  was  plotted  over  d.b.h.  and  tree  height  to  deter- 
mine what  samples  might  be  combined.    The  following  tabulation  indicates  the  initial 
species  and  National  Forest  combinations  used  to  develop  surface  area  equations,  and 
the  number  of  observations  in  each  combination.    The  plotted  data  also  indicated  that 
variance  about  a  regression  line  would  be  related  to  the  independent  variable,  DH. 

Species  and  National  Forest  Observations 


»  Douglas- fir  on  the  Lincoln  and  Tonto  National  Forests  107 

Mexican  white  pine  and  white  fir  on  the  Lincoln  National  Forest 

Blackjack  pine  on  the  Coconino  National  Forest  285 

,  Blackjack  pine  on  the  Lincoln  and  Tonto  National  Forests  395 

Yellow  pine  on  all  Forests  165 


Least  squares  regressions  of  surface  area  on  DH  were  fitted  for  the  initial 
species  groupings.     Observation  weights  proportional  to  1/(DH)2  were  used  to  obtain 
minimum  variance  estimators  of  regression  coefficients.     This  initial  estimate  was  based 
on  work  with  tree  volume  equations  (Cunia  1964),  which  found  that  the  square  residuals 
increased  linearly  with  the  square  of  the  independent  variable. 

Analysis  of  covariance  (Freese  1964)  was  used  to  test  for  differences  between  the 
surface  area  equations  for  blackjack  pine  on  the  Coconino  National  Forest  and  for 
blackjack  pine  on  the  Lincoln  and  Tonto  National  Forests.    This  analysis  (table  1)  in- 
dicated that  the  differences  were  not  significant  at  the  5  percent  confidence  level,  so 
the  data  from  the  three  Forests  were  combined  and  a  general  equation  for  blackjack  pine 
was  developed.    This  equation  was  next  tested  against  the  yellow  pine  equation  and  the 
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Table  1 . --Analysis  of  covariance:  comparison  of  coefficients  for  unforked  blackjack 
pine  equation  on  the  Coconino  National  Forest  and  unforked  blackjack  pine 
equation  on  the  Lincoln  and  Tonto  National  Forests 


Species  location 

Degrees  of  : 

Residual  : 

Residual 

:  Calculated 

and  items  tested 

freedom  : 

Sum  of  squares  : 
 1_  

mean  squares 

:  F 

Coconino 

283 

2,126.34 

7.5136 

-- 

Lincoln  and  Tonto 

393 

3,517.91 

8.9514 

Total 

676 

5,644.25 

8.3495 

Difference  for  testing 

slopes 

1 

2.94 

2.94 

0.3521  N.S, 

Pooled  values 

677 

5,647. 19 

8.3415 

Difference  for  testing 

levels 

1 

1.51 

1.51 

.1813  N.S. 

Single  equation 

678 

5,648.70 

8.3314 

I  able  2  .--AnaLys 

is  of  covariance:  comparison 

of  coefficients  for 

unforked  yellow  pine 

equation 

and  unforked  blackjack  pine  equation 

Subspecies  and 

Degrees  of  : 

Residual  : 

Residual  : 

Calculated 

item  tested 

freedom  : 

Sum  of  squares  : 

mean  squares  : 

F 

Yellow  pine 

163 

25,564.70 

156.8390 

Blackjack  pine 

678 

5,591.36 

8.2468 

Total 

841 

31 ,156.06 

37.0464 

Difference  for  testing 

slopes 

1 

695.44 

695.44 

18.772** 

Pooled  values 

842 

31,851.50 

37.8283 

difference  between  the  equations  was  significant  at  the  1  percent  confidence  level 
(table  2).     However,  an  equation  combining  yellow  and  blackjack  pine  was  also  developed 
for  users  of  these  equations  who  do  not  differentiate  between  yellow  and  blackjack  pine. 
This  must  be  regarded  as  an  equation  with  an  omitted  variable — growth  vigor.    As  Kmenta 
(1971)  demonstrated,  this  could  result  in  biased  estimates  of  the  regression  coeffi- 
cients.    This  situation  should  not  be  unduly  alarming.     The  bias  involved  is  a  bias  of 
the  estimated  coefficient  as  an  estimator  of  the  true  parameter  relating  surface  area  to 
DH.     So  long  as  the  equation  is  used  to  predict  surface  area  for  trees  which  have  the 
same  relationship  between  DH  and  the  omitted  variable,  predictions  will  be  unbiased, 
though  inefficient. 

Using  the  final  species-Forest  groups,  the  following  model  was  fitted  by  least 
squares  regression  to  further  examine  the  form  of  the  residuals: 

In  [(Actual  SA- Predicted  SA) 2  ]  =  cq  +  c1  In  (DH) 

The  coefficient  ci  in  this  equation  estimates  the  power  by  which  the  squared  residual 
increases  with  an  increase  in  DH.  The  calculated  values  are  presented  in  table  3. 
The  powers  of  DH  are  close  to  2  (the  power  which  was  assumed  earlier)  for  all  species 
except  yellow  pine,  which  is  closer  to  1. 
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Table  3. --Coefficients  for  equations  of  residuals  and  certain  related 

regression  statistics 


„  :  „  :        „       :     Standard  error  :      Mean  square 

Species  group  Cn  C\  _  _  ■,  fxicr>\ 

r  6      r  :  u        :  1     :  of  Ci  :  residual  (MSR) 


Ponderosa  pine  -10.78196  2.04002  0.092805  5.00786 

Yellow  pine  -4.56996  1.19232  .375083  5.13841 

Blackjack  pine  -9.75613  1.85458  .125818  5.04500 

Douglas-fir  -9.77705  1.83869  .220917  3.57956 

Mexican  white  pine,  white  fir           -9.98315  1.88056  .304233  6.45851 


An  estimate  of  the  mean  squared  residual  about  the  regression  surface,  as  con- 
ditioned by  DH,  can  be  obtained  from  the  equation: 


r  fI,m  MSRn 

S         =  e  tco  +  cl  <DH>  +  — 1 


yx 
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where    =  1/2  mean  square  residual,  which  is  necessary  for  conversion  of  logarith- 

2       mic  estimates  to  arithmetic  units  (Baskerville  1972). 

A  final  set  of  surface  area  regressions  was  then  fitted  with  observation  weights 
proportional  to  l/(DH)Cl.    These  equations  and  certain  regression  statistics  are 
presented  in  table  4.    The  equation  for  ponderosa  pine  was  forced  through  the  origin 
because  the  intercept  was  negative,  which  violated  the  expected  model  form.    All  other 
coefficients  meet  expectations. 

For  those  who  wish  to  calculate  estimated  confidence  intervals  for  predicted  values 
of  surface  area,  the  values  for  mean  weighted  DH  and  the  weighted  corrected  sum  of 
squares  have  been  included  in  table  4. 


Table  4 .--Surface  area  equations  and  certain  regression  statistics 


Surface 

area  : 

Sy -x 

Weighted 

:  Corrected 

Species  group 

equation  : 

(Weighted  root  mean 

mean 

:  weighted 

square  residual) 

DH 

:  sum  of  squares 

Ponderosa  pine 

0 

146150 (DH) 

3.2655 

242.064 

82,142,100 

Yellow  pine 

3 

10548  + 

149661 (DH) 

4.8301 

1,351.650 

83,183,300 

Blackjack  pine 

0 

52975  + 

143295 (DH) 

3.0411 

248.056 

21,302,200 

Douglas-fir 

3 

31473  +  . 

130655 (DH) 

3.8076 

360.927 

10,561,300 

Mexican  white  pine, 

2 

01727  +  . 

137988(DH) 

4.6509 

298.479 

9,370,900 

white  fir 
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ANALYSIS  OF  FORKED  TREE  DATA 


Forked  tree  data  were  grouped  by  species,  surface  areas  were  plotted  over  d.b.h.  and 
height,  and  resulting  trends  examined.     The  initial  grouping  of  data  and  the  number 
of  observations  for  each  group  were: 

Species  Observations 


Blackjack  pine  62 
Yellow  pine  15 
Douglas-fir,  Mexican  white  pine,  white  fir  19 

The  data  were  next  examined  by  calculating  the  ratio  of  the  actual  forked  tree  sur- 
face area  divided  by  the  predicted  surface  area  for  a  normally  formed  tree  and  plotting 
this  ratio  over  d.b.h.,  over  height  of  the  tallest  stem,  and  over  the  ratio  of  height  to 
fork  divided  by  height  of  the  tallest  stem.     No  relationship  was  apparent  over  d.b.h.  or 
over  height  of  the  tallest  stem,  but  there  was  a  relationship  over  the  ratio  of  height  to 
fork  divided  by  height  of  the  tallest  stem.     This  relationship  appeared  to  take  the  form: 

FSA  FH 
CPSA}  =  do  "  dl  (TH) 


where  FSA  =  forked  tree  surface  area,  in  square  feet 

PSA  =  predicted  tree  surface  area,  in  square  feet,  for  an  unforked  tree 
of  the  same  d.b.h.  and  height 
FH  =  height  to  fork,  in  feet 

TH  =  height  of  tallest  stems  in  the  forked  tree,  in  feet 


Physical  considerations  require  that  as  the  position  of  the  fork  approaches  the  tip  of 
the  tree,  the  predicted  forked  tree  surface  area  should  approach  the  prediction  for  an 
unforked  tree.     To  meet  this  requirement,  the  following  model  was  fitted  through  the 
origin  for  each  species  group: 


This  model  can  then  be  changed  back  to  the  original  form  through  these  transformations: 


d    =  1  +  f 

o 


di  =  (-f) 


Analysis  of  covariance  was  then  used  to  test  differences  in  coefficients  for  the 
three  species  groups.     No  significant  differences  were  found  at  the  5  percent  confidence 
level   (table  5),  so  the  data  were  combined  and  the  following  regression  was  fitted: 

-  1.214839  -  0.  214839  (^J) 


To  use  this  equation,  the  predicted  surface  area  for  a  normally  formed  tree  is 
calculated  using  the  equations  previously  presented  in  table  4.     This  predicted  surface 
area  is  then  multiplied  by  the  value  calculated  in  the  above  equation  to  estimate  the 
forked  tree's  surface  area.     For  those  users  who  do  not  measure  the  height  to  the  fork, 
the  following  equation  can  be  used: 

FSA  =  1.086481  (PSA) 

where  FSA  =  estimated  forked  tree  surface  area  in  square  feet.     This  equation 

FSA 

is  based  on  the  mean  value  of  C50T")  • 
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Table  5. --Analysis  of  covariance:  comparison  of  coefficients  for  forked  tree  equations 


Species  and 
items  tested 


Degree  of 
freedom 


Residual 
sum  of  squares 


Residual 
mean  squares 


Calculated 
F 


Blackjack  pine  60 

Yellow  pine  13 

Douglas-fir,  Mexican  white  pine,  17 

white  fir   

Total  90 

Difference  for  testing  slopes  2 

Pooled  values  92 


0.572876 
.301581 
.193569 


1.068026 

.032780 
1.100806 


0.0093914 
.0215415 
.0107539 


0.0118670 

.0163900 
.0119653 


1.3811  N.S. 
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